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ABSTRACT 

We measure the spatial clustering of metals in the intergalactic medium from z ~ 1.7 to 3.0, as traced 
by 643 C IV and 104 Si iv > 10^^ cm~^ absorption systems in 19 high signal-to-noise (40-80) and 
high resolution (i? = 45000) quasar spectra. The number densities and two-point correlation functions 
of both these species are largely constant with redshift, suggesting the bulk of metal ejection occurred 
at z > 3. However, at z < 1.9 some additional signature appears in the C iv correlation function 
at 500 km/s, associated with four strong and peculiar systems. At all redshifts, the C iv and Si iv 
correlation functions exhibit a steep rise at large separations and a flatter profile at small separations, 
with an elbow occurring at ~ 150 km/s. We show that these properties are consistent with metals 
confined within bubbles with a typical radius Rs about sources of mass > Mg, and use numerical 
simulations to derive best-fit values of i?^ 2 comoving Mpc and Mg ~ 5 x 10^^ at z = 3. This 
does not exclude that metals could have been produced at higher redshifts in smaller, but equally 
rare, objects. At the level of detection of this survey, IGM enrichment is likely to be incomplete and 
inhomogeneous, with a filling factor ~ 10%. 

Subject headings: intergalactic medium, quasars: absorption lines, cosmology: observations, large-scale 
structure of the universe 



1. INTRODUCTION 

No one knows where they came from, how they got 
there, or when it happened, but quasar (QSO) absorp- 
tion line studies have encountered heavy elements in all 
regions of the tenuous intergalactic medium (IGM) in 
which they were detectable (eg. Songalia & Cowie 1996). 
First, measurements of Nciw/Nyii indicated that typi- 
cally \C/H] ~ —2.5 in somewhat overdense regions at 
z ~ 3 , with an order of magnitude scatter (eg. Ranch, 
Haehnelt, & Steinmetz 1997). Later, statistical meth- 
ods pushed further, showing that C iv may be present 
in underdense regions (Ellison et al. 2000), and that a 
minimum IGM metallicity of approximately 3 x 1{)~^Zq 
is already in place at z « 5 (Songaila 2001; hereafter 
SOI). Although the presence of metals in underdense re- 
gions has not yet been firmly established (eg. Petitjean 
2001; Carswell et al. 2002; Schaye et al. 2003; Aracil et 
al. 2003a), their very existence outside galaxies has pro- 
found cosmological implications. 

As the long cooling times of large clouds of primordial 
composition are shortened by even modest levels of en- 
richment (Sutherland & Dopita; Scannapieco & Broad- 
hurst 2001), IGM metals accelerated the formation of 
massive (^^IO^^Mq) galaxies. Similarly, the violent 
events that ejected heavy elements from galaxies have im- 
portant implications for the thermal and velocity struc- 
ture of the IGM (eg. Gnedin & Ostriker 1997), and would 
have exerted strong feedback on dwarf galaxies forming 
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nearby (eg. Thacker, Scannapieco, & Davis 2002). 

Still, no one knows where they came from. Although 
numerous starburst-driven outflows have been observed 
at 3 < z < 5 (Pettini et al. 2001; Frye, Broadhurst, & 
Benitez 2002), it is unclear whether such objects are re- 
sponsible for the majority of IGM enrichment. In fact 
a variety of theoretical arguments suggest that these 
galaxies represent only the tail of a larger population 
of small "pre-galactic" starbursts that formed earlier 
(Madau, Ferrara, & Rees 2001; Scannapieco, Ferrara, & 
Madau 2002). On the other hand active galactic nuclei 
are observed to host massive outflows (Weyman 1997), 
whose contribution remains unknown. Finally, theoreti- 
cal studies have suggested that the primordial generation 
of metal-free stars were very massive (eg. Schneider et 
al. 2002), resulting in a profusion of tremendously pow- 
erful pair-production supernovae that distributed metals 
at z > 15 (Bromm et al. 2003). 

This paper describes the first direct attempt to use 
clustering measurements to discriminate between these 
disparate models. Making use of 643 C iv and 104 Si iv 
absorption systems, measured in 19 high signal-to-noise 
quasar spectra, we are able to place strong constraints on 
the spatial distribution of intergalactic metals at inter- 
mediate redshifts (1.5 < z < 3.1 and a detection limit of 
Nciv ^ 10^^ cm~^). Whatever the objects that enriched 
the IGM, it is clear that they formed in the densest re- 
gions of space, regions that were highly clustered. As this 
"geometrical biasing" is a systematic function of mass 
(eg. Kaiser 1984), the large-scale clustering of metal- line 
systems encodes valuable information on the masses of 
the sources from which they were ejected. Likewise, we 
are able to relate the small-scale clustering of metals to 
the typical radius enriched by each source. 

The structure of this work is as follows. In §2 we sum- 
marize the properties of our data set and reduction meth- 
ods. In §3 we present the number density of C iv and 
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Si IV absorption systems and their two point correlations. 
In §4 we compare these correlations with simulations to 
constrain the masses and ejection energies of the sources 
responsible for IGM metal enrichment. A discussion is 
given in §5. 

2. DATA SET AND REDUCTION 

2.1. Observations 

The ESO Large Programme "The Cosmic Evolution of 
the IGM" was devised to provide a homogeneous sample 
of QSO sight-lines suitable for studying the IGM from z 
= 1.7-4.5. High resolution {R ^ 45000), high signal-to- 
noise (40 and 80 per pixel at 3500 and 6000 A, respec- 
tively) spectra were taken over the wavelength ranges 
3100-5400 and 5450-9000 A, using the UVES spectro- 
graph on the Very Large Telescope. Although the com- 
plete emission redshift range is covered, emphasis is given 
to lower redshifts, and the C IV and Si IV metal lines dis- 
cussed in this paper were well-detected over the redshift 
ranges of 1.5-3.0 and 1.8-3.0 respectively, as described in 
Table 1. 

Observations were performed in service mode over two 
years. The data were reduced using the UVES context 
of the ESO MIDAS data reduction package, applying the 
optimal extraction method, and following the pipeline 
reduction step by step. The details of this reduction and 
spectrum normalization are described in further detail in 
Aracil et al. (2003a). 



Name 



C IV 



Coverage 

Si IV 



PKS2126-158 

Q0420-388 

HE0940-1050 

HE2347-4342 

HE0151-4326 

Q0002-422 

PKS0329-255 

Q0453-423 

HE1347-2457 

HE1158-1843 

Q0329-385 

HE2217-2818 

Q1122-1328 

Q0109-3518 

HEOOOl-2340 

PKS0237-23 

PKS1448-232 

Q0122-380 

HE1341-1020 



3.280 
3.117 
3.084 
2.871 
2.789 
2.767 
2.703 
2.658 
2.611 
2.449 
2.435 
2.414 
2.410 
2.404 
2.263 
2.222 
2.220 
2.190 
2.135 



2.36- 
2.23- 
2.21- 
2.04- 
1.97- 
1.96- 
1.91- 
1.87- 
1.83- 
1.71- 
1.70- 
1.68- 
1.68- 
1.67- 
1.56- 
1.53- 
1.53- 
1.50- 
1.46- 



3.24 2.74- 

3.08 2.59- 

3.04 2.56- 

2.83 2.38- 

2.75 2.31- 

2.73 2.29- 

2.66 2.23- 

2.62 2.19- 

2.57 2.15- 

2.41 2.01- 

2.40 2.00- 

2.37 1.98- 

2.37 1.98- 

2.36 1.97- 

2.22 1.84- 

2.18 1.81- 

2.18 1.81- 

2.15 1.78- 

2.10 1.74- 



3.24 
3.08 
3.04 
2.83 
2.75 
2.73 
2.66 
2.62 
2.57 
2.41 
2.40 
2.37 
2.37 
2.36 
2.22 
2.18 
2.18 
2.15 
2.10 



Table 1. List of lines of sight. Here Zcm is the quasar 
emission redshift. Only C iv and Si iv systems detected 
redward of the Lya forest and at least 3000 km/s blueward 
of Zcm are used. 

2.2. Line Identification 

Metal-line systems were identified using an automated 
two-step procedure. For each species we applied a pixel- 
by-pixel procedure which compared the spectrum with 
one that was rescaled and shifted according to the ratio of 
the oscillator strengths, extracting the minimal flux and 
accounting for blends. A detection threshold was then 
applied to these spectra, such that only absorption fea- 
tures with equivalent widths (EWs) larger than 5 times 



the noise rms were accepted. Finally a set of physical 
criteria were applied, to eliminate false detections. 

For each identified system the final decomposition into 
subcomponents was carried out using VPFIT (Carswell 
et al. 1987). For each component, initial z values were 
taken from the central wavelengths of the sub-features, 
while initial values for the Doppler parameter b and col- 
umn density N were computed directly from the EW and 
peak absorption of each sub-feature. This procedure is 
described in detail in Aracil et al. (2003b) and has been 
tested on simulated spectra, doing well for all systems 
with realistic values of N and b. 

Finally we applied a set of four cuts to the automated 
hst generated by VPFIT: log [cm^^j > 12 due to the 
detection limit of our procedure, 6 > 3 km/s to avoid 
false detections due to noise spikes, logiV[cm~^] < 16 
to remove two very badly saturated subcomponents, and 
6 < 45 km/s to avoid false systems due to errors in con- 
tinuum fitting. Our cuts resulted in a final data set of 
643 C IV and 104 Si iv systems, drawn from an original 
sample of 748 and 149 respectively. 

3. ANALYSIS 

3.1. Distribution Functions 

This sample was first used to compute the column 
density distribution function f{N). Following Tytler 
(1987), f{N) is defined as the number of absorbing sys- 
tems per unit column density and per unit redshift path 
X{z) = I [(1 -I- z)^/^ - 1] ■ In an ilm = 1 universe, which 
is a reasonable approximation over the redshift range 
considered, f{N) does not evolve for a population whose 
physical size and comoving space density are constant. 
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-111 



o 




-lyE 



12.0 12.5 13.0 13.5 14.0 14.5 15.0 

log N 

Fig. 1. — Column density distribution of C iv (upper panel) 
and Si iv (lower panel) absorption systems. In each panel, 
systems are divided into two redshift bins: 1.5 < z < 2.3 
(squares) and 2.3 < z < 3.1 (triangles). Column density 
bins are lO'^'^A'^ cm~'^ wide and ila error bars are given, 
computed from the variance of the data between spectra. 
The dashed lines are the power law fit measured in SOI. 

In Fig. 1 we plot f{N) for both C iv and Si iv sys- 
tems, dividing our sample into two redshift bins (1.5-2.3 
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and 2.3-3.1). Both species are consistent with a lack of 
redshift evolution, as found by previous studies of C iv 
and Si iv absorption systems (SOI; Ellison et al. 2000). 
Furthermore the overall density distribution of C IV is 
consistent with a power-law of the form f{N) oc N~^-^ 
as fit by SOI. From this figure, we see that some in- 
completeness sets in around 10^^'^ cm^^. While fewer 
in total, Si iv systems are also consistent with a lack of 
evolution, and follow a similar power law with a lower 
overall magnitude. 

3.2. Clustering of Metal Absorption Systems 

Having constructed a sample of well-identified metal 
absorption systems, whose distributions are consistent 
with previous measurements, we then computed their 
two-point correlation function in redshift space, S,{v). 
This quantity was last studied by Ranch et al. (1996) 
and Boksenberg, Sargent, & Ranch (2003) who noted a 
marked similarity between ^(u) of C iv and Mg II (see 
also Petitjean & Bergeron 1990, 1994). For each quasar, 
we computed a histogram of all velocity separations and 
divided by the number expected for a random distribu- 
tion. We then combined these distributions, weighting 
them according to the number of pairs in each quasar, 
yielding ^{v) + 1. The corresponding error bars were 
computed according to the appropriate weighted vari- 
ance between quasars. 
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Fig. 2. — Two point correlation function of C iv (upper 
panel) and Si iv (lower panel) absorption systems. In each 
panel our systems have been divided into two redshifts bins, 
with symbols as in Fig. 1. The solid line in the upper panel 
includes the four unusual low-redshift systems (see text). 

The resulting correlation functions are shown in Fig. 2, 
again split into two redshift bins. While the spatial distri- 
bution of both species shows no evolution at separations 
< 400 km/s, our comparison uncovers a low-redshift fea- 
ture in the 500 km/s C iv bin. This we have traced back 
to arising from four extended systems amongst our 19 
QSO sightlines, and seems to be due to some peculiar- 
ity in the internal structure of these systems. These will 
be studied in detail in a future publication, and are ex- 



cluded from our model comparisons below, although this 
cut has no effect on our final fits. 

The overall shape and amplitude of the C iv and Si iv 
correlation functions are strikingly similar and are con- 
sistent to within the Si iv measurement errors. Both 
functions exhibit a steep rise at large separations and 
a flatter profile at small separations, with an elbow oc- 
curring at 150 km/s. Unlike Lya absorption systems 
(see Cristiani et al. 1997), the correlation of C iv does 
not seem to depend strongly on absorption column den- 
sities. While choosing subsets below a fixed N increases 
the noise of our measurements, the shape of ^{v) is pre- 
served to within error bars. 

4. CONSTRAINTS ON METAL SOURCES 

The two-slope shape of the measured correlation func- 
tions is highly suggestive of a picture in which metal 
bubbles of a typical size (associated with the small scale 
clustering) are generated about objects of a typical mass 
(whose geometrical bias is associated with the large-scale 
clustering) (see also Petitjean & Bergeron 1990). In or- 
der to explore this connection further, we generated a 
simple model, based on a dark-matter-only simulation. 
As no evolution was detected, and the properties of both 
species were similar, we focused our attention on com- 
paring the properties of C iv absorbers with simulated 
groups at a single representative redshift of three. 

We adopted a Cold Dark Matter cosmological model 
with h = 0.7, no = 0.3, = 0.7, Ob = 0.05, cTg = 0.87, 
and n = 1, where ilg, ^^A, and Of, are the total matter, 
vacuum, and baryonic densities in units of the critical 
density (pc), cr| is the variance of linear fluctuations on 
the 8/i~ Mpc scale, and n is the "tilt" of the primor- 
dial power spectrum. We made use of a simulation with 
512^ dark matter particles, contained with a box 50/ h 
comoving Mpc on a side. The mass of each particle was 
7.7 X 10^ Mq, which gives a minimum mass resolution for 
our group flnding of 10^'^ Mq as we select only groups 
with 40 or more particles. The initial conditions were 
produced using GRAFIC (Bertschinger 2001) and were 
used as the inputs to the parallel version of GADGET 
(Springel et al. 2001). Halo detection was performed us- 
ing the HOP algorithm (Eisenstein & Hut 1998) with the 

cuts (5pcak = 160, ^saddle = 140, and (5outor = 80. 

The resulting list of groups was then used to construct 
a simple simulated picture of IGM metal enrichment. We 
selected all 2; = 3 groups above a threshold group mass 
Ms as markers of the centers of enrichment events, about 
which we painted spheres of radius Rg. Note that our 
choice of this set is not meant to imply that enrichment 
occurred at z = 3, but rather that it occurred at an un- 
known redshift higher than the observed range, centered 
on groups whose large-scale clustering was equivalent to 
z — 3 objects of mass Mg. All smaller groups within such 
spheres were associated with C IV absorbers, each with 
an overall impact parameter equal to a factor b times its 
dark matter radius. The observed number of C IV sys- 
tems per unit redshift path was well-reproduced by fixing 
b = 10. 

Thirty-two planes were taken, and for each plane all 
512 sightlines were drawn; the average two-point cor- 
relation was then compared with detected ^ values, as 
shown in Fig. 3. As expected, the large scale correla- 
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Fig. 3. — Comparison between detected C IV correlation functions (points) and simulated correlations (lines). In each panel Ms is fixed, 
Rs is varied over three representative values, and b = 10. 



tions are almost completely dependent on Ms while the 
position of the elbow in ^(u) is set by Rg. Our mea- 
surements are clearly suggestive of i?s — 1.6 Mpc and 
M = 5 X IO^^Mq z = 3 objects, corresponding to a fill- 
ing factor of / ~ 10%, at the level of detection of our 
survey. Varying b over a wide range of values has only 
a minor impact on this comparison, changing the inner 
slope of ^(v) slightly while giving a much poorer fit for 
the overall number of systems. 

A straightforward geometric argument allows us to as- 
sociate Rs with a wind ejection energy per object, Es. 
Here we adopt a simple picture in which the average bub- 
ble expands according to a Sedov- Taylor solution (Sedov 
1959). In this case Eg is related to the bubble radius 
as Eg = 48tt R^ pt~^ / 125 , where t is the expansion time 
and R is in physical units such that R = i?s/(l + Zo), 
where is the redshifts of observation, respectively. If we 
assume for simplicity that t is approximately the Hubble 
time (2.3 Gyrs at 2 = 3 in our cosmological model) and 
that the shell expands into material at the mean density 
at this redshift, this gives E = 1.6 x lO'"^^ erg, which is 
equal to the kinetic energy input from roughly 1.3 x 10^ 
Typell supernovae (SNIIe). This is roughly consistent 
with the overall number of C IV systems, which can be 
converted into a cosmological density of fie « 5 x 10"^ 
as described in SOI. 

5. DISCUSSION 
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We have studied the clustering of 643 C iv and 104 
Si IV systems detected in UVES observations of 19 high-z 
QSOs using an automatic procedure. The independently 
derived correlation functions of C iv and Si iv are similar, 
giving confidence in our automated approach. The shape 
of the correlation function is consistent with a picture 
where metals are confined within bubbles of i?s ~ 2 Mpc 
about Ms ^ 5 X 10^^ M© halos at z = 3. This imphes 
that the filling factor of the metals is only 10% at the 
detection limit of the survey. 

As our overall mass constraint is derived only from 
the geometrical bias of the sources, however, we can- 
not exclude models in which equally rare, but somewhat 
smaller higher-redshift objects are responsible for metal 
enrichment. Furthermore, the observed number density 
of C IV shows little evolution at 2 < 3 as at higher red- 
shift (eg. SOI). While more detailed investigations are 
necessary, it remains that all such models will be strongly 
constrained by the clustering of intergalactic metals. 
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